/* this creates the models used to create Table 1 and Table 2 in the main text */

import delimited "WHC_replication_british_only.csv"

global X   indirect_rule groupsize lnarea lncarea lnpop lngdp indep 
global Z0  coastdist 
global Z0a coastdist_indirect

* Model 1
probit status_egip $Z0 $X if sample1==1, cluster(cowid) nolog
* Model 2
probit status_egip $Z0 $Z0a $X if sample1==1, cluster(cowid) nolog

biprobit (status_egip= $Z0 $Z0a $X) (family_onset_fe_flag= status_egip $Z0 $X) if sample1==1, cluster(cowid) nolog